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Abstract 

QCD and QED exhibit an infinite set of tliree-point Green's functions that contain 
only OZI rule violating contributions, and (for QCD) are subleading in the large 
expansion. The Green's functions describe the "decay" of a J^^ = {1,3,5...} ^ 
exotic hybrid meson current to two J = (hybrid) meson currents with identical P 
and C. We prove that the QCD amplitude for a neutral hybrid {1, 3, 5 . . .} ^ exotic 
current to create r/vr'' only comes from OZI rule violating contributions under certain 
conditions, and is subleading in A'^c- 
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1 Introduction 

More than a decade ago Lipkin argued that an explicitly identified contribution to the decay 
of a J^*" = {1, 3, 5 . . .}^^ "exotic" hybrid meson to rjn vanishes |l[]. Here J denotes the in- 
ternal angular momentum, P (parity) the refiection through the origin and C (charge conju- 
gation) particle-antiparticle exchange. These are conserved quantum numbers of the strong 
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and electromagnetic interactions. By "hybrid meson" we mean a fermion-antifermion with 
additional gauge bosons. "Exotic" means that the J^^ cannot be constructed for conven- 
tional mesons in the quark model, or equivalently that there are no local currents built only 
from a fermion and antifermion field that can have these J^"-^. 

Lipkin's intuitive argument was later extended and placed in a more formal context, called 
"symmetrization selection rules" 0. However, these need to be converted into rigorous 
quantum field theoretic arguments, which is the subject of the current work. 

We outline Lipkin's intuitive argument for the decay of a positively charged = 1~, 3~, . . . 
state to TjTi^ when strong interactions with isospin symmetry is assumed. G-parity con- 
servation implies that the neutral isospin partner of the initial state is J^^ exotic. We 
consider the decay process where a quark and antiquark in the initial state proceed in such 
a way that the quark ends up in the one final meson, and the antiquark in the other meson 
(Fig. 111). This decay process is called "connected", because each meson is connected to 
the other mesons via quark lines. The gluons are not indicated. The fact that the neutral 
isospin partner of the initial state is J^*-^ exotic, and that the initial state contains a quark 
and antiquark, implies that the initial state is a hybrid meson. The argument is depicted in 
Fig. m. Taking the initial hybrid at rest, the rj and vr"*" emerge with momenta — k and k re- 
spectively. First consider the three top left diagrams. The top diagram has a negative sign 
in front by convention. When the transformation k ^ — k is applied, the middle diagram 
is obtained, noting that the decay is in an odd partial wave, which acquires a minus sign 
under the transformation. This is a general property of odd partial waves. The bottom 
diagram is obtained by noting that the amplitude to create a uu pair is the same as for a 
dd pair by the assumption of isospin symmetry, which treats the up and down quarks the 
same. The three top right "hadronic" diagrams are now obtained from the three top left 
"quark" diagrams by attaching the initial hybrid to the initial ud quarks, and the final vr"*" 
to the final ud quarks. Since the flavour wave function of the r] is proportional to uu + dd, 
it is attached to either uu or dd, with a positive relative sign. Because each of the three top 
left quark diagrams are equal, it follows that each of the three top right hadronic diagrams 
are equal. The two bottom diagrams depict the decay amplitude, taking into account that 
there are two possible ways for the final rj and tt^ to couple, since the quark in the initial 
state can go either to the rj or the vr. Looking back at the top right hadronic diagrams one 
immediately notices that the decay amplitude vanishes. This is the desired result. 
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Figure 1: Lipkin's intuitive argument. 

The argument above can be repeated without assuming isospin symmetry, as will be the 
case in the remainder of this work, if the initial state is neutral. Hence isospin is not an 
essential assumption. It appears strange that a decay amplitude vanishes from very general 
considerations if the decay is allowed by the conserved quantum numbers of the strong 
interaction. This paradox is resolved when one notices that it was only argued that the 
connected contribution to the decay vanishes, not the entire decay amplitude. Lipkin's 
argument serves as a guide to how a quantum field theoretic argument would proceed. 

In this Paper we demonstrate that some explicitly identified contributions to certain three- 
point Green's functions vanish. This is called field symmetrization selection rules. Par- 
ticularly, the connected (Okubo-Zweig-Iizuka (OZI) allowed) contribution to each Green's 
function vanishes. Hence, the Green's function only has a disconnected (OZI forbidden) 
contribution, which is expected to be phenomenologically suppressed by virtue of the OZI 
rule The Green's function is built from two J = (hybrid) meson currents with identical 
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P and C, and a J^*" = {1,3, 5...} ^ exotic hybrid meson current. 

We subsequently investigate the physical consequences. The amplitude for a {1, 3, 5 .. .}""*" 
hybrid current to create t^tt" is shown to be proportional to the Green's function under 
certain conditions. Because some explicitly identified contributions to the Green's function 
vanish, it follows that the amplitude does not come from these contributions, called sym- 
metrization selection rules at the hadronic level. Particularly, the amplitude does not arise 
from the connected contribution to the Green's function, and hence only from the discon- 
nected contribution, which is expected to be suppressed by virtue of the OZl rule. The 
experimental consequences pertain to a central issue in hadron spectroscopy: the search for 
hybrid meson bound states beyond conventional mesons and baryons. 

For Quantum Chromodynamics (QCD) with a large number of colours N^., the foregoing 
conclusions will be expressed as follows. Because a disconnected contribution to a Green's 
function is subleading in the large Nc expansion, and the Green's function only has a discon- 
nected contribution, the Green's function is subleading in Nc- Also, since the amplitude for 
a {1, 3, 5 . . .}~^ hybrid current to create r/vr*^ only comes from the disconnected contribution 
to the Green's function, it is subleading in A*",,. 

It has previously been shown that the connected part of the quenched Euclidean three-point 
Green's function of specific hybrid meson neutral 1~"^ and two pseudoscalar ( J = 0) currents 
vanishes exactly in QCD, if isospin symmetry is assumed We remove the three italized 
superfluous assumptions, but still need the connected part and hybrid meson currents. We 
also generalize beyond specific currents, beyond 1""*" and pseudoscalar currents, beyond the 
connected part, and beyond QCD. 

In section 2 the currents and Green's functions are introduced. The principle of symmetriza- 
tion is developed. In section 3 the Green's functions are calculated leading to succinctly 
stated field symmetrization selection rules. An explicit example is discussed. In section 4 
the physical consequences are investigated, yielding symmetrization selection rules at the 
hadronic level, which are concisely stated. Section 5 studies the selection rules in the large 
Nc expansion. Section 6 contains further remarks. 
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2 Symmetrization 



Consider local currents of the form 



{Af,{x),B{x),C(x)} = ^tpiix) {aiaf,{x),bib{x),Cic{x)} ipiix) 



(1) 



where ipi{x) is a quark or lepton field of flavour i, and ai,bi,Ci are c-numbers weighting 
the flavours. The currents are diagonal in flavour. The "matrices" a^(a;), b{x) and c{x) 
contain an arbitrary number of Dirac matrices, Gell-Mann colour matrices, derivatives 
(acting both to the left and the right), gluon or photon fields and correspond to gauge 
invariant currents A common choice for the flavour structure of the three currents is 
uu — dd, uu + dd — ss and uu — dd, interpolating for an isovector resonance, an rj and a 
7r°. The matrix a^{x) in Eq. |l] is chosen to ensure that P = —. We require that currents 
B and C have J = 0, as well as equality of b{x) and c{x), implying equal P and C. Note 
that equality of the matrices b{x) and c{x) does not imply equality of the currents B{x) 
and C{x). Because C is the same for both these currents, charge conjugation conservation 
requires that the C of A^^^z) is +. Take to have odd J 0. Hence A^ must be chosen 
to have J^*" = {1, 3, 5 . . .} ^. These are exotic quantum numbers, and is built from a 
fermion and an antifermion field, so that the current is a hybrid current, i.e. has to contain 
at least one gluon or photon field. For the current with J 7^ 0, the appropriate Lorentz 
indices are indicated by fi. 

We start by demonstrating that certain three-point Green's functions are equal to their 
antisymmetric parts. This is done by first defining the spatial Green's function, G^{x, y, z), 
and decomposing it into symmetric and antisymmetric parts. The Green's function of 
interest will be the Fourier transform of the spatial Green's function, G^(p,t). We then 
argue that the Fourier transform of the symmetric part of the spatial Green's function, 
G^{p,t), vanishes. Hence, G^(p,t) equals the Fourier transform of the antisymmetric part 
of the spatial Green's function, G^(p,t). 

Define the three-point Minkowski space Green's function 



with xq = Ho = t, Zq = and z = 0. The Green's function describes the decay (production) 



G,{x,y,z) = {0\B{x) C{y) A,{z) |0) 



(2) 
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process of a current Afj^{z) at time propagating into the final (initial) currents B{x) and 
C{y) at some positive (negative) time t. Although we shall refer to G'^( "Green's 
function", the usual usage of the term requires the currents to be at different times, i.e. 
Xq ^ Ho, with the currents ordered from positive to negative times. The Green's function 
Gf^{x, y, z) can be written as G^^{x, y, z) + Gf^{x^ y, z), i.e. the sum of parts symmetric and 
antisymmetric under exchange of x and y, as any function of x and y can be written. 

Define the Fourier transforms 

G^(p,t) = I d^xd'ye'^^-^-^-y^G^ix,y,z) (3) 

{G^(p,t), G;^(p,t)} = J d^xd'ye'^^^-^-^y^ {G'^{x,y,z),G^{x,y,z)} (4) 

Note that these Fourier transforms cannot be inverted to give the spatial Green's functions, 
since they only have one momentum variable p. From Eqs. ^ and ^ it follows that ^^(p, t) = 

G^(p,t) + G;^(p,t). 

By exchanging integration variables x ^ y 

G,i-p, t) = J d'x d'y e^(P-P-y) G^((y, t), (x, t), z) 
= J ci^^ d^y e^(P™) {G'^{y,x,z) + G^{y,x,z)} = G'^{p,t)-G^{p,t) (5) 

Exchanging integration variables x — x and y — y yields 

G,i-p,t) = -G,ip,t) (6) 

by using that G^((— x, t), (— y,t), z) = —G^{x,y, z), by conservation of parity. We hence- 
forth restrict ourselves to the parity conserving theories of Quantum Electrodynamics 
(QED) and QCD. We also used that the product of the parities of the currents B{x), 
C{y) and A^{z) is —1, which follows from the assumptions below Eq. |l[ This means that 
the decay (production) process is in odd partial wave. Eq. ^ is a well-known property of 
such a process. 

Combining Eqs. ^ and ^ yields the result 

Gl{p.t) = Q (7) 
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Bose symmetry states that identical bosons are not allowed in an odd partial wave. The 
way one shows this would be analogous to the steps above if G^{p,t) = 0. Hence, the 
Fourier transform of the symmetric part of G^{x, y, z) in Eq. |^ vanishes by arguments that 
are the field theoretical version of Bose symmetry. 

From Eq. |^ follows the desired result 

G>,t) = G;^(p,t) (8) 



3 Field Symmetrization Selection Rules 

We proceed to explicitly identify some contributions to ^^(p, t) that vanish. This is attained 
via partially evaluating G^(p,t). It is subsequently shown that the action of a certain 
operator, Op, on the Fourier transform of the antisymmetric part of the Green's function. 
Op G^(p,t), does not contain some contributions. Because G^(p,t) = G^(p,t) from Eq. 
H, it follows that Op G^(p,t) does not contain these contributions. 

We now evaluate the contributions to the Green's function G^(x, y, z) with the currents in 
Eq. H, using that Gf^{x, y,z) = \ y, z) - G^{y, x, z)). 



G^^{x,y,z) = -{0\ {B{x)C{y)-B{y)C{x)) A,{z)\0) 

= ^ H^iCj (0| ( ^i{x) b{x) ^i(x) i!j{y) c{y) ^j{y) - 

^,{y) h{y) ^,{y) ^,{x) c{x) ^,{x) ) A^{z) |0) (9) 

The various contributions to this expression are now discussed. Consider contributions to 
the expression where the same flavours are isolated in the currents B{x) and C{y), i.e. 
contributions where i = j. Using Eqs. |^ and P, these contributions to G^{p,t) can be 
written 



i 

X (o| [ h{x) , Uy) Ky) i^^{y) ] M^) |o) = Kip,t) (lo) 



where we used b{x) = c{x), and denoted the contributions by A^. The important ob- 
servation is that the difference of currents in Eq. ^ has simphfied to the commutator of 
currents in Eq. It is possible to show that has a polynomial dependence on p 

(see Appendix) 0. There hence exists a polynomial operator (containing a derivative of 
high enough power in p) with the property that Op A^(p,t) = 0. This result will later be 
demonstated in an explicit example. We conclude from Eq. |10|that Op G^{p,t), and from 
Eq. 1^ also Op G^(p, t), do not contain contributions from the same flavour in currents B{x) 
and C{y), a result which corrects the former treatment [§]. Although this is the desired 
result, and the end of the mathematical derivation, it will be pivotal to develop a more 
intuitive understanding of the contributions. 

Op G^{p,t) is an operator acting on the Fourier transform of G^{x,y, z). For the purpose 
of illustration in the next few paragraphs, consider the time t to be positive, with B{x) 
slightly advanced at time xq = t + 5t with respect to C (y) at time yo = t in the definition of 
Gfj,{x, y, z) in Eq. 0. Here 5t is small, positive and non-zero. The sign of t or the magnitude 
of 5t will not change our eventual conclusions. Since all the currents are at different times, 
with the currents ordered from large to small times, it follows that G^{x,y,z) is a Green's 
function according to the usual usage of the term, which implies that it can be represented 
by a path integral. For concreteness, consider the contribution to G^{x, y, z) in Eq. ^from 
the up quark flavours in the currents B{x), C{y) and v4^(z) in Eq. |I], which is (modulo 



VA Vip Vip 5{f{A)) deiM-F u{x)b{x)u{x) u{y)c{y)u{y) u{z)a^{z)u{z) exp(— i j d^x C) 

= j <^(/(^)) detA^iT det7o5-^ exp(-z j d^x C^^^) 
{TT[-a,{z) S^{z,x) b{x) S^{x,y) c{y) S^{y,z) - a,{z) S^{z,y) c{y) S^{y,x) b{x) S^{x,z)] 
+ Tr[a^(^) S^{z,x) b{x) S^{x,z)] Tr[c(i/) St{y,y)] + Tr[a^(z) St{z,y) c{y) St{y,z)] 
X Tr[6(x) S^{x,x)] + Tr[a^(^) S^{z,z)] Tr[6(a;) St{x,y) c{y) St{y,x)] 
- Tr[a,{z) S^{z,z)] Tr[6(a:) S^{x,x)] Tr[c(y) S^{y,y)]} (11) 

The gauge-fixing condition 6{f{A)) and Faddeev-Popov determinant detj^p are indicated. 
The QCD / QED Lagrangian is denoted by C and the part containing terms without 
fermions C^'^^ The up quark propagator in a background field is S:^{x, y) and the fermion 
determinant (containing fermion loops) det7o5'~^ 0. 
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(a) (b) 
Topology 1 Topology 2 Topology 3 

Figure 2: Topologies contributing to Eq. pT| . 



The terms on the right hand side (R.H.S.) of Eq. [TT| correspond to the topologies in Fig. ^. 
This is seen by associating with each up quark propagator an up quark line in the figure. 
Also, z is associated with the left hand side (L.H.S.) of each topology, and x, y with the 



blobs on the R.H.S. In Eq. |Tl| the first two terms correspond to the quark line "connected" 
topology 1 in Fig. ^, the second two to topology 2, term five to topology 3b and term six 
to topology 3a. 

The topologies hence represent all the different ways that a fermion and antifermion field 



on the L.H.S. of Eq. |TT] can be "contracted" by the fermion integration to yield the fermion 
propagators on the R.H.S. 

Similar manipulations to Eq. |ll| can be performed for contributions where not all the 
fermions in the currents are up quarks, or not all the fermions have the same flavour. 

From Fig. ^ it can be seen that topologies 1 and 3b have the interesting property that 
when a certain flavour contributes from current B{x), then the same flavour contributes 
from current C{y). This is because the currents are diagonal in flavour. Although B{x) and 
C{y) contain various different flavour structures according to Eq. |l|, topologies 1 and 3b 
force only the same fermion flavours in B{x) and C{y) to contract. For the other topologies, 
it is sometimes the case that the flavours from two currents are the same, but not always. 
Recalling that Op G^{p, t) does not contain contributions from the same flavour in currents 
B{x) and C{y), we derive 

Field symmetrization selection rules (FSSR): All contributions to the three-point 
Greenes function Op G'^(p,t) from the connected topology 1 and topology 3b (and some 
from topologies 2 and 3a), i.e. contributions from the same flavour in currents B{x) and 
C{y), vanish exactly for all momenta p and times t for an infinite set of equal matrices 
b{x) = c{x) and different flavour structure currents B{x) and C{y) in QCD and QED. The 
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Green's function decribes the "decay" of a J^^ = {1, 3, 5, . . .}~'^ exotic hybrid meson current 
to two J = (hybrid) meson currents B{x) and C{y), which are identical {b{x) = c{x)), 
e.g. have identical parity and charge conjugation, except possibly for their flavour. 

The contributions that vanish by the FSSR are those from the same flavour in currents 
B{x) and C{y). Since the product B{x) C{y) occurs in the Green's function (see Eq. |^), 
these are contributions of the form biCi ipi{x) b{x) ipi{x) 4'i{y) c{y) ipi{y), where i indicates 
the flavour, and no summation is implied. It is evident that the precise values of 6j and q, 
i.e. the flavour structure of B{x) and C{y), do not affect the fact that the contributions 
vanish by the FSSR. The same is true for v4^(z). For example, both the common choice 
of uu — dd, uu + dd — ss and uu — dd for the flavour structure of Afj,{z), B{x) and C{y), 
and the alternative choice of respectively uu — 2dd, uu + dd — 2ss and 2uu — dd, obey the 
FSSR. 

The FSSR do not depend on the various parameters of QCD and QED, i.e. masses, cou- 
plings, charges, number of colours and flavours, but do require the CP violating 6 parameter 
to be zero, since this parameter violates parity conservation. FSSR also occur is pure QED, 
where QCD interactions are turned off. 

We now discuss an example to which the FSSR apply. The gauge invariant isovector- 
like local current A^{z) = {T'^j,)^ (u{z)F^j^{z)X''-fiu{z) - d{z)F!^j^{z)X''jid{z)'^ is J^^' = 1"+ 
exotic 0, where 7j and A" are Dirac and Gell-Mann colour matrices respectively, and 
{T^jk)^ is a tensor combining the spatial indices i,i,k to build a spin 1 object. It is a 
hybrid current, since it contains a gluon field, as can be seen by the presence of the gluon 
field tensor F°'j^. The gauge-invariant isoscalar- and isovector-like local currents B{x) = 
'u(x)75'u(x) + J(x)75c/(x) — s(x)75s(x) andC(?/) = u{y)'y5u{y) — d{y)'y^d{y) are pseudoscalar 
{J^'-'' = 0^^), with 75 a Dirac matrix. Note that b{x) = c{x) = 75, but that B{x) 7^ C{x). 
Also, we do not assume isospin symmetry. To evaluate in Eq. |10|, the commutator 
[u{x)'y5u{x) , u{y)-y5u{y)] must be evaluated. The commutator can be shown to equal 5^(x — 
y) ( m(x)757o75u(?/) — ^i(l/)757o75^(a;) ). Because of the delta function, it is immediate that 



A^(p,t) in Eq. |T0| is independent of p, i.e. it has a polynomial dependence on p. The 
polynomial operator Op = has the promised property that Op A^(p,t) = 0. In this 
example, the commutator actually vanishes because the delta function forces x = y, so that 
5^(x - y) ( m(x)757o75m(i/) - M(l/)757o75M(a;) ) = 0; implying that A^ = 0, and one can 
take Op = 1. There does, however, exist examples where Op is not trivial as in this case. 
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For example, one can show that replacing b{x) = 75 by b{x) = ''i^F^^i^x) F"-^^^{x) yields a 
non-trivial Op. 



4 Symmetrization Selection Rules at the Hadronic Level 

We now investigate the physical consequences of the FSSR, particularly the amplitude for 
a {1,3,5...}^^ hybrid current to create r^vr^, which can be obtained from the Green's 
function. This is done via an alternative route to the former treatment Q], which contains 
erroneous aspects [[10|] . 

A natural quantity to obtain from the three-point Green's function is the amplitude for 
the hybrid current to create a stable two-body state. The hybrid current is not expected 
to interpolate for a stable particle, so that we shall not be able to extract the T-matrix 
for a stable hybrid particle to decay to the two-body state. To obtain the amplitude, a 
complete set of asymptotic, i.e. stable, states n are inserted in the Green's function. The 
leads to a general relation (Eq. |14D between the Green's function and the amplitudes for a 
{1, 3, 5 . . .}^^ hybrid current to create the asymptotic states. Under certain conditions this 
equation can be simplified to show that the Green's function is proportional to the amplitude 
for a {1, 3, 5 . . .} ^ hybrid current to create t^tt^, Eq. [1^. Because some explicitly identified 
contributions to Op 0^(p, t) vanish, it follows that the amplitude does not come from these 
contributions. 



Restricting to QCD on its own |]rT|, from the definition Eq. 



G^(p,t) = / rf^xrf^^e^^P-^-Py)^: (0|i?(x,0) O(y,0) e^*^» (n|A^(^)|0) (12) 

where we used time translational invariance of the fields B{x) and C{y), e.g. i?(x, t) = 
e^^*5(x, 0)e-^^*, with H the QCD Hamiltonian. The product fi(x,0) O(y,0) and 
should be colourless to make the expression non-zero, which is the case since each current 
has been assumed to be gauge invariant. If the quarks have their physical masses, the 
lowest asymptotic states contain the states vr and 1], the lowest stable states of the QCD 
spectrum . The vr only decays weakly and electromagnetically, and is hence stable as far 
as QCD on its own is concerned. The t] has a width of only 10^^ times the typical hadronic 
width of 100 MeV, so that it is very nearly stable. 

Space translational invariance of the fields B{x) and C{y), e.g. i?(x, 0) = e~*^'^i?(0)e*^ ''. 
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with P the QCD momentum operator, is now employed. Performing one of the integrations 
in Eq. [12| (the one over x + y) 



^^(p,^) = (27r)3 53(pj e-'E^t ^Q|( / g.p.x ^(^^Q)) ^(Q)|^^ {n\A^iz)\0) (13) 

where the integration variable x — y is denoted by x. Here p„ and are the momentum 
and energy of state n. The Euclidean space analogue of all the steps up to here is obtained 
by taking t —it. However, we now restrict to Minkowski space in order to integrate over 
time. From Eq. [l^ 



/ rft6pG^(p,t)e*^* 

J —oo 

= (27r)^ 5\vn) 5{En - E) Op (0|( / d\ e^^'^ i?(x, 0)) C(0)|n) (n|A^(^)|0) (14) 

where is a real number. The delta functions indicate that the asymptotic states are at 
rest and have energy E. 

For E below the tt+tt^tt^ or tt^tt^tt^ (called the Trvrvr) threshold only vrvr asymptotic states 
contribute. Since for two-pion states at rest (vr^k 7r° — k|A^(z)|0) and (vr+k vr^ — k|y4^(2;)|0) 
respectively vanish by Bose symmetry and CP conservation, this forces the L.H.S. of Eq. 



14| to be zero. If E is above the titttt threshold the sum in Eq. |1J can be shown to be 



infinite. However, we shall show below that if the r/vr" threshold can be made below the 
TTTTTT threshold, in contrast to experiment, the sum has only one contribution. This is 
important because we want a single amplitude on the R.H.S. of Eq. |1^ to be proportional 
to the L.H.S. 

In QCD we are free to tune the quark masses away from the masses corresponding to 
experiment. Taking the r^vr" threshold to be below the vrTrvr threshold for some range of 
quark masses, is called the "QCD dynamics" condition. This condition can heuristically 
be satisfied by noting that four times higher up / down quark masses would yield a two 
times heavier vr, while the rj mass would not change much, so that rj can move below the vrvr 



threshold. This follows from the fact that the vr and rj masses, and m^, are c^viiu + vn^ 



and cy (4ms + + mrf)/3 by chiral symmetry breaking, with mu,d,s the current quark 
masses and c a constant. Under the QCD dynamics condition t] lies below tttt threshold, so 
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that it is stable. This means that the formahsm is exact when asymptotic states involving 



T] are inserted in Eq. 12 



Consider E between the rjn^ and tttttt thresholds. Then the sum in Eq. ^is over all on-shell 
TT and 1] states with momenta ki and k2 respectively, i.e. X)n = / d^ki /(27r)^ / d^k2 /(27r)^. 



Noting that pn = ki+k2 and En = ykf + m^ + y'kj + m^, and performing the integrations, 
Eq. O implies that 



JEt ^ 



-oo 



dt Op G„(p, t) e^^* = f{q') V(^' - + ^v)')iE' - (m. - m,y 



X (E^ - (m^ + m^)2(m^ - m^)^) J dn^, Qm Op {0\i J d^x e'^ '' 5(x, 0)) C(0) |7rk r/ - k) (15) 

Here q'^ = ki — k2 and ki = — k2 = k. In the expression is a function of |k|, and |k| is 
restricted by En = E, so that is a function of E. The amplitude for a hybrid current to 
create r] and vr mesons, {nk r] — k\A^{z)\0) , should occur in Eq. |15|. However, we used the 
Lorentz properties of (vrk r] — k\Ani{z)\0) to define it as qm.fi'f), with m a spatial index, 
for having J = 1 |13|. For a time index m and other J similar results follow. 

On the R.H.S. of Eq. |r^, (0|-B(x, 0) C(0)|7rk ?7 — k) is a property of r] and vr. For example, 
it contains a term proportional to (0|i?(0)|7rk) (0|C(0)|?7k), a product of individual prop- 
erties of TT and ?7, when the currents -B(O) (C(0)) have the same quantum numbers as the 
asymptotic states tt (r/). One hence interprets Eq. |T5| as meaning that the decay described 
by the three-point Green's function on the L.H.S., is proportional to the decay amplitude 
described by f{q^) on the R.H.S. , up to a "constant" of proportionality which describes 
properties of r] and vr. Suppression of the decay on the L.H.S. should hence translate into 
suppression of the decay on the R.H.S., since the properties of r] and vr should remain 
unaltered. 

From Eq. |T5|, noting that f{q^) is proportional to the L.H.S., 

Symmetrization selection rules at the hadronic level (SSR): The amplitude for 
a neutral J^^ = {1,3,5...} ^ exotic hybrid meson current to create or annihilate rjTT^ , 
Qmf{q^), does not arise, in QCD, from contributions to Op G.m(p,t) that vanish by the 
FSSR. This holds for quark masses chosen such that the rjir^ threshold is below the tt+tt^tt^ 
or 71^71^71^ thresholds, and for E between these thresholds. 

This concludes the physical consequences of the FSSR. 
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5 Large Nc 



We now study the effect of taking A^^^ to be large in QCD . Here one makes a classification 
based on power counting in Nc- The amplitude (yrk 77 — k|A^(2)|0) is 0{1) |15|, so that f\q^) 
is 0{1). The contributions to the Green's function from topology 1 is 0{Nc) to leading 
order, from topologies 2 and 3b 0{1) and from topology 3a 0{^) so that the L.H.S. 
of Eq. ^ would ordinarily be 0{Nc)- However, since topology 1 does not contribute by 
the FSSR, the L.H.S. is and hence f{q^) is 0{^) ||16[. Hence both quantities are 

subleading to their usual Nc counting. Whence 

^6pG'^(p,t)^0 /(g')--O asiV.^oo (16) 

where both expressions ordinarily remain non-zero as A^^^ 00. The two exactly vanishing 
expressions are respectively the consequences of the FSSR and SSR in large Nc. The SSR 
in large Nc is: The hybrid J^'" = {1,3,5. . .}^^ ?77r° amplitude is 0{^) in QCD, i.e. 
vanishes exactly in the large Nc limit, for quark masses chosen such that the rjir^ threshold 
is below the tt'^h'tt^ or tt^tt^h^ thresholds, and for E between these thresholds. This is not 
the same as Bose symmetry, since r] and 7r° are not identical particles in the large Nc limit. 
The finding that the amplitude for the hybrid current to decay to 777?° is subleading in Nc 
intuitively follows from the fact that OZI-violating processes are large Nc suppressed. 

The preceding discussion in this section assumed that the L.H.S. of Eq. [1^ is "ordinarily" of 
0{Nc). This is attained when (0|5(x, 0) C(0)|7rk - k) on the R.H.S. of Eq. |l| is 0{Nc). 
This is true when the currents B{x) {C{y)) individually have the same quantum numbers as 
the asymptotic states vr {rf) or rj (tt) ||T5|, i.e. when the currents are pseudoscalar. Up to this 
point in the Paper we have not required either of the currents to have a specific parity and 
charge conjugation. However, if the currents are not pseudoscalar, (0|-B(x, 0) C(0)|7rk r] — k) 
is at most 0{\/Wc) [l^, so that the L.H.S. cannot be 0{Nc), and hence should be regarded 



as ordinarily 0{1): the next possibility in the Nc counting. The FSSR still yield that 
the L.H.S. is 0(1), but as this is not new information, the FSSR imply no additional 
restrictions on f{q^). Thus the SSR in large Nc are only interesting when B{x) and C{y) 
are pseudoscalar currents. 

This concludes the statements of the FSSR and SSR in large Nc- A few final remarks are 
in order. 
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6 Remarks 



Firstly, on the relationship between SU{3) and large N^. For equal mass up, down and 
strange quarks {SU{3) flavour symmetry) the t] and vr are among the degenerate lightest 
states of QCD, satisfying the QCD dynamics condition that the r/vr^ threshold is below 



the 711171 threshold, implying the existence of SSR. The SSR and Eq. |15| yield that the 
hybrid J^*^ = {1,3,5. . .} ^ i]7c^ amplitude vanishes exactly because topologies 2 and 
3a in Fig. ^vanish in the SU{3) limit, noting that an SU{3) octet current B{x) or C{y), 
interpolating respectively for t] or 7r°, does not couple to a quark-antiquark pair created 
from the vacuum. It has been known independently for some time that the SU{3) octet 
{1,3,5...}^+ rj7i^ amplitude vanishes exactly in the SU{3) limit 0, |T^. 

Hence the hybrid {1,3,5. . .}""'" r^vr^ amplitude vanishes exactly in either the large N^. 
or SU{3) limits, but the one does not follow from the other, as SU{3) symmetry does not 



derive from, or does not imply, the large Nc limit |T^. The amplitude should be more 
suppressed than either limit indicates, due to the constraints from the other limit. 

Secondly, on the role of Bose symmetry. One finds that contributions to the Green's func- 
tion that vanish by a field theoretical version of Bose symmetry, at least after the polyno- 
mial operator is applied, do not contribute to the hybrid {1, 3, 5 . . .}^+ — ^ ?77r° amplitude. 
However, this amplitude does not itself vanish by Bose symmetry. 

We now remark on the field symmetrization selection rules. The vanishing contributions 
to the Green's function were foreshadowed by, and have direct analogues in, the "sym- 
metrization selection rules I" of the non- field theoretic analysis of ref. 0, where decays 
of {1,3,5...} hybrids to two J = (hybrid) meson states which are identical in all 
respects except possibly flavour are prohibited. The latter condition is translated into the 
requirement that h{x) = c{x) for the FSSR. In ref. the selection rule applied to flavour 
components of the (hybrid) mesons B and C which are identical, e.g. uu for both. This is 
exactly the case for the contributions to the Green's function for which we have FSSR. 

The symmetrization selection rules at the hadronic level have important consequences for 
models. The hybrid {1, 3, 5 . . .}~^ i]7t^ amplitude only comes from disconnected topolo- 
gies. This feature puts the SSR in contradiction with most current models of QCD, al- 
though most find it in an approximate form. Particularly, the flux-tube and constituent 
gluon models all find an approximate selection rule for the connected decay of the low-lying 
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1^+ hybrid to rjir^ [0. In these models the decay is proportional to the difference of the 
sizes of the i] and vr wave functions, which is in contradiction with the SSR. Constituent 
gluon models also have low-lying hybrids called "quark excited" hybrids whose connected 
decay to rjir^ vanishes exactly, consistent with the SSR [Q. A non-zero decay 1 ^ rjn^ 
via final state interactions has been estimated from the decay of 1 ^ to two mesons which 
then rescatter via meson exchange to r/vr^ [|19[. The process is described in QCD by con- 
nected decay (with a quark loop), so that it contradicts the SSR. In practical calculations 
in the above models the QCD dynamics condition may not be satisfied, so that the model 
calculations are strictly not required to obey the SSR. However, model parameters can be 
changed to satisfy the QCD dynamics condition, so that the SSR have to be obeyed. Hence 
models that do not incorporate vanishing connected topologies are inadequate. In QCD 
sum rules, the connected topology vanishes, consistent with the FSSR; and /(g^) is small, 
consistent with expectations from the OZI rule [Q. The QCD sum rule calculations |20| 
constitute explicit examples of the results of this work. 

If one's aim is to obtain information about the physically interesting hybrid {1, 3, 5 . . .} ^ 
TjTT^ amplitude, it is possible to do so from a variety of Green's functions. The essential 
point about a quantum field theory approach is that information about amplitudes can 
be extracted from various Green's functions. In addition to Green's functions involving a 
fermion-antifermion current going to two fermion-antifermion currents considered in this 
Paper and elsewhere Green's functions containing a quark-antiquark current going 

to a pure glue and a quark-antiquark current have been considered in QCD sum rules pO| . 
The latter case can be shown to yield no selection rules for Green's functions (FSSR), and 
hence cannot be used to deduce selection rules for physical amplitudes (SSR). However, the 
fact that the SSR cannot be deduced does not imply that the SSR is not valid. 

We developed SSR for the hybrid {1,3,5...}^+ rjix^ amplitude. Related SSR have 
been found for four-quark and glueball initial states 0. These SSR are derived in a 
formal context which can transparently be extended to a rigorous quantum field theoretic 
argument along the lines of this work. Extension to two-body final states beyond r/vr^ is 
more conceptually involved, due to the analogue of the QCD dynamics condition. In the 
interest of brevity these extensions should be considered in future work. 

Lastly, we remark on experimental consequences. The Crystal Barrel experiment has re- 
cently claimed evidence for pp —>■ resonant 1^^ rjir^ at a level that is a significant 
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fraction of especially the P-wave pp annihilation [^. Although the branching fraction of 
1 ^ — s> rjTT^ is not known, it is reasonable to assume that it is substantial. This is qualita- 
tively at odds with the SSR if the l""*" resonance is interpreted as a hybrid. Even though 
this is not rigorous, as the QCD dynamics condition is not satisfied experimentally, the 
small change in quark masses from their experimental values needed to enable the validity 
of the QCD dynamics condition indicates that the 1 resonance is qualitatively incon- 
sistent with being a hybrid meson. Other possibilities for the interpretation of the 1 
enhancement have recently been discussed in refs. [Ol 111, O, |23 . 
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Appendix 



We show that A^ defined in Eq. |TU| has a polynomial dependence on p if the number of 



derivatives acting on individual fields in the currents B{x) and C{y) is bound. 

The commutator in Eq. [10| can be expressed, by a general property of commutators, as a 
sum of terms, each of which contains either a commutator of two boson (gluon or photon) 
fields, or an anticommutator of fermion ip and conjugated fermion ip (quark or lepton) 
fields. Here we used the fact that boson fields commute with fermion fields, that two 
fermion fields anticommute, and that two conjugated fermion fields anticommute. Each 
of the (anti) commutators is proportional to delta functions 5^(x — y) (or derivatives act- 
ing on delta functions), by virtue of the canonical (anti) commutation relations of fields 
at equal time t. For example, for commutators of photon fields [y4^(x, t), 74j^(y, t)] = 0, 
[y4^(x, t), Aj^ly, t)] = igfj_y5^{^ — y), . . .; and for anticommutators of lepton and conjugated 
lepton fields {^^^(x, t), ^<x(y, t)} = -flJ^i^ - y), {V'c(x, t),i>^{y, t)} = -7 ■ d^5^{x - y), • • 
Spacial derivatives and dy might be acting on these (anti) commutators, since b{x) in 
Eq. |l^ will in general contain derivatives acting on fields. A derivative dy can be expressed 
in terms of a derivative 9x since the delta function only depends on x — y. The possi- 
bility of temporal derivatives acting on the fields have already been incorporated in the 
(anti) commutation relations. Hence a generic term contributing to Eq. |10| is of the form 



I d'x d'y e^(P— P-y) /^(x, y, z) d^5\^ - y) (17) 

One now performs integration by parts over the variable x, which yields powers of p when 
the derivatives act on the exponential, as well as derivatives acting on fn{x,y,z). (There 
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is no surface term as the delta function does not contribute for x far from y). Eventually 
there will be no derivatives acting on the delta function. When performing one of the 
integrations, the delta function forces e**-P ''~P'^^ = 1, so that the only p dependence is the 
various powers of p. Hence has a polynomial dependence on p. This is true as long as 
the number of derivatives n in Eq. |1^ is bound. 
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